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Abstract 

We compute in the Ginzburg-Landau approximation the gluon Meissner masses for the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) phase of QCD with three flavors in the kinematical range where 
it is energetically favored. We find real Meissner masses and therefore chromomagnetic stability. 
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I. INTRODUCTION 



The attractive interaction between quarks in the color antisymmetric channel leads at high 
densities and small temperatures to their Cooper pairing, see 0] and for reviews U 0] ; 
in particular, for three flavors at asymptotically high densities, the Color-Flavor-Locking 
(CFL) phase is the ground state of the theory, characterized by a spinless color- and flavor- 
antisymmetric diquark condensate 0. 

The physically more interesting pre-asymptotic phases, where the mass of the strange 
quark and the chemical potential differences S/i due to /? equilibrium cannot be neglected, 
are subject of intense study at this moment. Proposals include the 2SC phase j2|, the gapless 
phases g2SC jfj, and gCFL 0,0. Unfortunately the appearance of imaginary gluon Meissner 

makes the gapless phases instable, and also the 
possible antidotes to cure the chromo-magnetic 



masses (for g2SC see [9|, for gCFL see 
2SC phase shows instability 
instability are discussed). 



gCFL see 



However, for appropriate values of Sji, quarks may form pairs with non-vanishing total 
momentum: p x + p 2 = 2q ^ 0, see K| and for a review l^j ]. which leads to a Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) 15 1 phase. For two flavors it has been shown Q| that 
the 2SC instability implies that the LOFF phase is energetically favored; however, it is not 
clear if the neutral LOFF phase can cure the chromomagnetic instability of the two flavor 
superconductive quark matter [l^, [3| • 

LOFF with three flavors is however the much more difficult but physically interesting 
case at intermediate densities. A first study of the problem was carried out in a Ginzburg- 
Landau (GL) approximation It was found that condensation of the pairs u — s and 

d — u is possible in the form of the inhomogeneous LOFF pairing. 

The problem which comes next is then to find out whether such a phase is chromomagnetic 
stable. This is the subject of the present note. Within the GL expansion we find that this 
is indeed the case. This result is another indication that the LOFF phase of QCD plays 
a major role at intermediate hadronic densities. Since these pre-asymptotic densities are 
probably relevant for the cores of compact stars, a whole field of investigation is opened up, 
with items ranging from transport properties, cooling processes and the glitches in the pulsar 
rotational frequency to the possible implications for gravitational radiation during collisions 
of compact objects and the mechanisms of gamma ray bursts. Besides, given the universality 
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which is peculiar to the phenomena occurring at the Fermi surface, one expects important 
indications to come from the study of the similar problems with cold fermionic atoms 2o| . 
mainly because of the greater parameter flexibility possible in laboratory experiments. 

The plan of the paper is as follows. In section |H] we discuss the model and review the 
results of Q]; in section HTT1 we present the computation of the gluon Meissner masses; in 
section HVl we discuss our results and in section IS we give our conclusions. The concluding 
appendix contains some useful formulas. 



II. THE THREE FLAVOR LOFF PHASE OF QCD 

In the CFL phase of QCD all the eight gluons acquire real Meissner masses, while in- 
troducing the strange quark mass and the neutrality constraints (which leads to the gCFL 
phase) some of these masses become imaginary, which is a signal of instability. In this paper 
we compute gluon Meissner masses in the three flavor inhomogeneous case (LOFF phase). 
We will work in the GL approximation For a system of two massless (u and d) and one 
massive (s) quarks the QCD action is: 
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d*x (i ptf - M°f + $ 7o ) • (1) 



Here a, /3, are color indexes, i, j flavor indexes; M°f = S a/3 diag(0, 0, M s ) is the quark mass 



matrix and Df- = d x 5 a/3 5ij + igA a T^Sij; ji^a is the matrix of the chemical potentials, see 
below. It depends on \x (the average quark chemical potential), /i e (the electron chemical 
potential), and /13, fi$, related to color [71 \21\. 



We treat the strange quark mass at its leading order, ie by a shift in the strange quark 
chemical potential u s — > u s — and we adopt the High Density Effective Theory 

(HDET) |22j (see |^ and I23J for reviews). HDET uses the fact that, for T -> and at 
weak coupling, the relevant modes of the QCD action are those near the Fermi surface. It is 
useful therefore to decompose the quark momentum as follows: p = fin + £; fin is the large 
component (n a unit vector representing the quark Fermi velocity) and £ is the small residual 
momentum that one can take in the direction of n using reparameterization invariance, so 
that £ = £n. Moreover one introduces velocity dependent fields ip n and \l/ n corresponding 
to positive and negative energy solutions of the Dirac equation: 

1>{*) = J |p MnX bM*) + *„(*)) ' ( 2 ) 



Substituting (J2J) in (JTJ) and integrating out the negative energy components, one gets in 
momentum space, at the next-leading-order in 





a(3\ 


y>e + 2n_ 


ij / 



(3) 



where (P)?? = - gA a Tfb ih with I* = (p°, fn). Moreover V = (l,n), V = 

(1, -n), /Zjf = (j$ - US'* Sij - M s 2 /(2/i) 5 a %5 i3 and = - (v»V v + 1>V™) /2 . 

In the color-flavor and Nambu-Gorkov (NG) basis introduced in JJJ] the free propagator 
reads 



S 



[S?]ab 
[S ( 22 



JAB 



(Po - £ + A 4 a) 1 

Oo + ^-AtA) -1 



(4) 



where /i a = (firu,fi g d, fibs, fird,figu, firs, fibu, figs, fibd)- Writing the propagator as in Eq. (jU) 
one doubles the fermion modes, which is compensated by an extra factor 1/2 in momentum 
integration. 

To keep into account quark condensation we add to the QCD action the bilinear quark 
term 

1a = -\ J d'x j ^A°f(x) ^,- n (x)C 75 ^,„(x) + h.c. (5) 
where the gap function is given by [19| 



A ? = E A ' exp < 2zqi ■ x > ea/3/ e ^ • 



(6) 



i=i 



Eq. © corresponds to assume a Fulde-Ferrell ansatz for each inhomogeneous pairing; 2qi 
represents the momentum of the Cooper pair. The inverse fermion propagator in the super- 
conductive phase becomes 



o-l 
D 



V 



AB 



A AB = -^2A I Tr[Fje I F B € I \ . 



(7) 



I>>o 2 ]ab j 



i=i 



We assume in this paper that global color neutrality is reached with almost vanishing 
color chemical potentials, 113 ~ /is ~ 0. This approximation is justified because the phase 
transition from the superconductive to the normal state is second order (see below), and the 
color chemical potentials are suppressed by inverse powers of /i. Therefore 
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(8) 



where Q is the quark electric-charge matrix. 

This formalism was applied in where a GL expansion of the free energy Q was 
performed. We briefly review here the results of this paper. From the ansatz (jHJ) it is clear 
that Q depends on 10 parameters, ie the three gaps, three q's and the electron chemical 
potential. The energetically favored state must be a global minimum of Q in the space of 
A's and q's, and has to be electrically neutral. As a consequence, one has to solve 

dtt dtt dn 

Wr ' a5T°- (9) 

The electrical neutrality condition leads to the result /i e ~ M 2 / (4/i) , which corresponds 
to a symmetric splitting of the Fermi surfaces of d and s around the u Fermi sphere. As a 
consequence Sfidu = 8(j, us = Sfi and Sfj,^ = 25fi. Moreover, for the computation of the gap 
parameters, it is sufficient to consider the condition dQ/dqi = at the 0(A 2 ), which leads 
to the well known relation q = 1.1997|5/i|. This implies that A2 = A3 and q2 = q3, which 
can be interpreted as the consequence of a symmetric splitting of the Fermi surfaces, and 
Ai = 0, due to the larger mismatch among the d and s Fermi surfaces. These results hold 
in the range (128 — 150) MeV for the parameter M s 2 //i 19j . For smaller values the LOFF 
free energy is higher than the gCFL energy; for higher values the normal phase is favored. 
For this solution the GL density of free energy reads 

n^ + S^ + ft^ + l^, do) 

where Q n is the usual normal contribution to the free energy; the explicit formulas for the 
coefficients can be found in the appendix. The transition to the normal phase at M 2 /fi pa 150 
MeV is second order, which justifies the GL expansion. 



III. MEISSNER MASSES IN GINZBURG-LANDAU APPROXIMATION 

In the effective theory there are two vertices describing the coupling of gluons and quarks. 
Either one gluon couples to a quark and an antiquark (three-body vertex, coupling ~ g) or 
two quarks couple to two gluons (four-body vertex, coupling ~ g 2 ). The two vertices come 
from terms in Eq. Q with one or two momenta £. At the order of g 2 the four-body coupling 
gives rise to the contribution g 2 n 2 /(2ir 2 ) identical for all the eight gluons; this result for the 
LOFF phase is identical to those of the normal or the CFL case (this term is independent 
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of the gap parameters). The three-body coupling gives rise to the polarization tensor: 

ilC(z,y) = -Tr[iS(x,y)iH^iS(y,x)iH v h ] (11) 

where the trace is over all the internal indexes; S(x, y) is the quark propagator, and H£ 
is the vertex matrix in the HDET formalism which can be read from Eq. (J3J). The quark 
propagator S has NG components (i, j = 1, 2), see Eq. (JJJ). At the fourth order in A: 

S 11 = S n + S n A [S 2 2 A* (S n + S n A S 22 A* S 11 )} , (12) 
S 21 = S 22 A* (So 11 + S 11 A S 22 A* S n ) , (13) 

where Sq 7 can be read in Eq. (@J); S 12 and S* 22 are obtained by the changes 11 «-> 22 and 
A <-> A*. 

|= a^t^^ vva= "^■ X/W CT) WN 

J = /W^)w\ = W^^WA = -^-XAAT ^^ VNA 

L= /w^^wv 

N= AAT^^WN = AAf^^^ 

FIG. 1: Diagrams appearing in the GL expansion of the gluon self-energy in Eq. Wavy 
lines denote gluons, solid lines are either Sq 1 and 5q 2 . Full and empty circles (squares) denote 
respectively the insertion of A2 (A3) and A2 (A3). Each A* is preceded by Sq 2 and followed by 
Sq 1 . For a detailed description see the text. 

The different contributions from the GL expansion are proportional to A , A 2 , A 4 . The 
contribution independent of A is equal to — g 2 /i 2 /(27r 2 ) and therefore cancels out the term 
from the four-body vertex. The 0(A 2 ), 0{A i ) terms arise from classes of diagrams with 
different topologies. Some of them are depicted in FigCJ The remaining diagrams are 
obtained first by duplicating all the diagrams by the exchanges (Sq 1 , A) <-> (Sq 2 , A*), then 



duplicating the diagrams presenting both the A2 and the A3 insertions by the exchanges 
A2 «-> A3, and the diagrams with unequal number of insertions on the quark lines by the 
exchange of the upper and lower lines of the fermion loop. 

The Meissner masses are tensors with spatial components and a nontrivial color structure: 

(x 2 )>-n> = o,p = o), (14) 

where i,j (a,b) are spatial (adjoint color) indices. We find 

M 2 im = 2 I l3 A\ + 2 J y A| + A3A2 (3 Kij + M l3 - 2 Ay , (15) 

M% m = 2 J y (A* + A3) + 2 J^Aj + A 4 3 ) + A3A2 (K^ + 2 M l3 ) , (16) 

M\ 33 = 2 I i:j A 2 2 + 3 JyA* + A3A2 (3 + M 4j - 2 Ny) , (17) 

^•,38 = -^/<,-A* -V3JijAi + AlAl( ^K ij -^M i j-^Nij + ±=L i ^j , 

(18) 

^J l88 = hj Qa 2 + + Jy (A} + 4A*) + A* A* ^ + I My - I A>y - | (19) 

(the expressions of /, J, K, L, M and N, represented in Fig. ^ are given in the ap- 
pendix); moreover the following identities hold: M 2 j22 = > 77 = j 
Mlj M = M 2 ijM (A 2 A 3 ) , A4?- )55 = -Mf i44 . Since the mass matrix of the gluons 3 
and 8 is not diagonal in the indices of the adjoint representation, one introduces the mass 
eigenstates A a = cos #^3 + sin 9iA iS , A i8 = — sin 9iA i3 + cos^A^; the physical masses are 
the corresponding eigenvalues. 

IV. RESULTS 

A. The limit cases A2 = or A3 = 

The limit A2 = is equivalent to the two flavor LOFF phase considered in Q] , with 

{ipai C75 ifrpj) oc A 3 exp {2zq 3 • r} e a p 3 e ij3 . (20) 
7 



In this limit we reproduce the results of confirming that for small gap parameters 
the Meissner masses of the screened gluons are real. First of all, for A2 = one has 
■M-ij a = M%j 22 = = 0' which is expected because there is an unbroken SU(2) color 

subgroup. Moreover M 2 ijM = M 2 j55 = M 2 jfi6 = M 2 j77 , where (A 3 = A): 

o 2 u 2 A 4 o 2 u 2 A 2 

xxM 96tt 2 (ql - 5/i 2 ) 2 ' ZM 8vr 2 (g c 2 - 5^) ' 1 J 

and q c = 1.1997 5fi given by a'(q c ) = 0. Finally, as in Eqs. (93) and (94) of [171] . 

In the same way one can treat the limit A 3 = 0: there are 3 massless and 5 massive gluons 
(with the same masses as before). The massless particles correspond to the generators T 4 , 
T 5 , f = T3/2 + V3T 8 /2 of the unbroken SU{2) subgroup. 



B. The three flavor LOFF case 



In the three flavor case, for A2 = A3, we obtain the results shown in Fig. |21 For the 
numerical evaluation we have used the results of Ref. [3] for the gap A 2 = A 3 . 




130 135 140 145 150 130 135 140 145 150 



Mg/H (MeV) mJ/h (MeV) 



FIG. 2: On the left: longitudinal squared Meissner masses, in units of the CFL Meissner mass, vs 
M 2 1 fj,\ from top to bottom the lines refer to the gluons A3, A§, A$ (long dashed), and A\ (dotted 
line). On the right: transverse squared Meissner masses; from top to bottom the lines refer to the 
gluons Aq, Ai, A3, and A&. 



On the left panel in Fig. El we report the longitudinal (ie zz) components of the squared 
Meissner masses against M 2 j\x, in units of the CFL squared mass J2J,|25( at the C(A 4 ) for 
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the representative gluons; on the right panel the results for the transverse (ie xx) squared 
Meissner masses are given. In both cases we obtain positive squared Meissner masses for all 
the gluons. The masses not reported are obtained from those displayed, according to the 
discussion at the end of section II VI 

From the figure it is clear that for each value of the strange quark mass, the transverse 
mass of a gluon is smaller than the longitudinal one. This is so because the transverse mass 
is zero at the 0(A 2 ) (see the appendix). As a consequence, it gets contribution only from 
the 0(A 4 ) and is therefore suppressed as A 2 /Sfi 2 in comparison with the longitudinal one. 
This behavior is analogous to the two flavor LOFF phase considered in Moreover, we 
find that the transverse mass of A s , although positive, is almost zero, being three order of 
magnitude smaller than the other ones. 

We conclude therefore that the LOFF phase with three flavors in the Ginzburg-Landau 
limit has no chromomagnetic instability. 

V. DISCUSSION AND CONCLUSIONS 

Some final comments are in order. First of all, we have computed the Meissner masses 
only for the single plane wave Fulde-Ferrell (FF) structure, see Eq. (JSJ). However from the 
two flavor case we know that more complicated crystalline structures have a lower free energy 
than the FF state and we expect the same to be true in the three flavor case j^ j. In the 
general case one should replace (JHJ) with 



where q} (i = 1,...,N) are the momenta which define the LOFF crystal; the geometry 
of the structure and the number N of plane waves should be determined by minimization 
of the free energy. Once the optimal structure is found, one should compute the Meissner 
masses. If this structure contains at least three linearly independent momenta, the Meissner 
tensor should be positive definite for small values of A, since it is additive with respect to 
different terms of (jHJ) to order A 2 fl^ |. These considerations suggest that a LOFF crystal 
can remove the chromo-magnetic instability of the homogeneous superconductive phases of 
QCD, resulting as the true vacuum of the theory. A second comment is that the results 
should be extended beyond the GL expansion. Finally, in a recent paper it has been 



3 




(23) 



i=l 1=1 
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found that, at strong coupling, gapless phases can be magnetic stable, and this direction 
also needs to be explored. 

In conclusion we have computed the gluon Meissner masses in the three flavor LOFF phase 
of QCD using the High Density Effective Theory and the Ginzburg-Landau approximation. 
The use of the GL expansion is justified in the LOFF window, 128 MeV < M^f/j, < 150 
MeV ^| . In this region we find that all the squared gluon Meissner masses are positive and 
therefore the LOFF phase of three flavor QCD is free from chromomagnetic instability. 



APPENDIX A: DEFINITION OF INTEGRALS 

In section imi we have introduced the following integrals: 

/..- f ^n n- I (Al) 

lJ ~ 4vr3 J An^^J (po-e + ^) 3 -(Po + ^-/i M + 2q-n) ' [Ai) 

,#V f dn f dp d£ 



4tt 3 J 4tt 3 J (p - £ + fin) 4 ■ (po + £ - fis + 2q • n) • (p + f - /2«j + 2q ■ n) ' 

(A3) 



T -^V [fa [ dpoMj ( h.A\ 

~ l 4^ J Att^^J ( Po + e-/i u -2q.ii)2-( P o-e + /id) 2 -(Po-e + /i-) 2 ' 



n ~ 4vr3 y 4vr ninj y ( P o-e+^) 3 -(Po+e-^ + 2q-np.( Po -e+/i,) ' lAbj 



AT.. - _,^! / ^-n n- ! ( M ) 

V ~ 4vr3 y 4vr n4 ^y ( Po + e-/i u + 2q-n)3.( Po -e + /i (i ) 2 -(po-e + ^) ' 1 J 

depicted in Fig. ^ To evaluate these integrals first we perform a Wick rotation po = ip 4 . 
The integral in £ is computed by a contour integration; to evaluate the integral in p 4 we 
note that the £ integration introduces the sign factor sign(j> 4 ) and the integrands depend 
only on ip 4 , therefore one can use 

/oo />oo 
dp 4 sign(p 4 ) F(ip4) = 2i^sm / <ip 4 F(zp 4 ) ; (A7) 
■oo </o+ 
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the remaining angular integrals are performed trivially. The prescription embodied in Eq. 
iTj) follows from the limit T — ► once one passes from finite to zero temperature. 
The integral iy = Iij(q, Sfi) has to be handled with care. In the case of I xx one has 



with a given by [3] 
a(q, 5[i) 



g 2 1 da(q, S/j) 
64 q dq 



(A8) 



4/x 2 



7T 
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■f 


— - log 


4(g 2 - 


) 
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- 5fi 


2 5 







(A9) 



the condition dQ/dq = at 0(A 2 ) reads o/(<7) — 0: thus vanishes at 0(A 2 ). Defining 
q c the value of q satisfying the equation a'(q c ) = and expanding around q c one finds 

,2 



64 q c 



-a"(qc) ; 



(A10) 



using dfl/dq = at the second order in A 2 one gets 



A 



<x"(q c )(q - q c ) = — (/?'(&) + WM) 



(All) 



where |12| (3 = fi 2 /[ir 2 (q 2 - 5fi 2 )}, and 
P23- 



(A12) 



7 47r (2q 3 • n + ^ u - fi d - ie) (2q 2 ■ n + ji u - fi s - ie) 

In (jAll|) the derivative of P23 is in the variable ^3, for fixed q2\ at the end one puts qi = q%. 
At the lowest non- vanishing order in A 2 one has 

<? 2 A 2 



64 2q c 



(/3'(? c )+2/4(g c )) 



(A13) 



The integral in the longitudinal treated in a similar way, 

I zz (q) « LMc) + (q~ QcKM ■ (A14) 
As a consequence, the integrals 1^ give contributions at the second and fourth order in A. 
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